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Abstract. We discuss the extension of radial SLE to multiply 
connected planar domains. First, we extend Loewner's theory of 
slit mappings to multiply connected domains by establishing the 
radial Komatu-Loewner equation, and show that a simple curve 
from the boundary to the bulk is encoded by a motion on moduli 
space and a motion on the boundary of the domain. Then, we show 
that the vector-field describing the motion of the moduli is Lips- 
chitz. We explain why this implies that "consistent," conformally 
invariant random simple curves are described by multidimensional 
diffusions, where one component is a motion on the boundary, and 
the other component is a motion on moduli space. We argue what 
the exact form of this diffusion is (up to a single real parameter k) 
in order to model boundaries of percolation clusters. Finally, we 
show that this moduli diffusion leads to random non-self-crossing 
curves satisfying the locality property if and only if k — 6. 



1. Introduction 

In this paper we discuss conformally invariant growing random com- 
pact sets in multiply connected domains. Our results are meant to 
provide some of the steps to extend the stochastic Loewner evolution 
of Schramm, ^Hj, from simply connected domains to multiply con- 
nected domains and Riemann surfaces. Based on Loewner's theory 
of slit mappings, Schramm showed that conformally invariant mea- 
sures describing random "simple" curves in simply connected domains, 
can be encoded into a diffusion on the boundary of the domain if the 
random simple curves satisfy a "consistency condition." Furthermore, 
under an additional but very natural symmetry condition, he showed 
that the diffusion on the boundary is a multiple of Brownian motion if 
the random simple curve lives in certain standard domains. 

We show that consistent and conformally invariant random simple 
curves in multiply connected domains can be encoded into a multidi- 
mensional diffusion. One component of this diffusion corresponds to 
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a motion on the boundary of the domain and the other components 
are the moduli of the domain with the random simple curve, grown 
up to time t, removed. The random simple curves we consider in this 
paper connect the boundary to the bulk (interior) and so the appro- 
priate moduli space is the moduli space of n-connected domains with 
one marked interior point and one marked boundary point. We show 
that consistency (a Markovian-type property) and conformal invari- 
ance essentially determine the diffusion up to the drift of the motion 
on the boundary. We call such diffusions on moduli space Schiffer dif- 
fusions. Under a symmetry condition, familiar from percolation, this 
drift component can also be identified, leaving a single real parameter 
k. Beginning with the Schiffer diffusion with this drift we show that the 
resulting family of random growing compact sets satisfies the locality 
property if k = 6. 

The fundamental observation that diffusion processes on the moduli 
space of bordered Riemann surfaces with marked points, given its path- 
wise solutions agree with the geometric constraints, yield the most gen- 
eral way to define probability measures on (simple) curves on surfaces 
and therefore contain "ordinary" SLE as a special case, was introduced 
in |H] . However, the current article is the first constructive imple- 
mentation of it, for the radial case and multiply connected domains. 
For general Riemann surfaces and the chordal case one proceeds along 
similar lines as described here, but with some necessary modifications. 
An important role then is played by the so-called "Hilbert Uniformisa- 
tion". 

The paper is structured as follows. In Section |H] we introduce a suit- 
able family of standard domains and describe canonical mappings onto 
these domains in terms of the Green function and associated functions. 
In Section |3] we introduce an appropriate time parameter for the Jor- 
dan arcs in a multiply connected domain, namely the conformal radius, 
and establish a variational formula for "increments" of the conformal 
radius under perturbations of the domain. In Section we establish 
what we call the radial Komatu-Loewner equation, which generalizes 
the radial Loewner equation to multiply connected domains. In Section 
El we obtain the corresponding equation for the evolution of the moduli 
and then show that the vector field in the differential equation satisfies 
a Lipschitz property, Theorem 16.11 This result is crucial as it allows 
us to reverse the construction: Start with a motion on the boundary, 
solve the equation for the moduli for the given boundary motion, then 
solve the radial Komatu-Loewner equation to obtain a growing family 
of compacts. In Section [7| we use the correspondence between growing 
"simple" curves and paths on moduli space to show that consistent 
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conformally invariant random simple curves are given by diffusions on 
moduli space, and identify the diffusion up to a single drift term. Un- 
der an additional symmetry condition we then identify the diffusion up 
to a single parameter k. In Section |H] we study the growing random 
compact sets for this diffusion and show that they satisfy the locality 
property if and only if k = 6. 

As this paper neared its completion we became aware of the thesis 
of Dapeng Zhan, which contains another version of the radial Komatu- 
Loewner equation, and discusses a class of random Loewner chains on 
Riemann surfaces, [T9] . The results we present here were previously 
announced in pQ. 

The authors would like to acknowledge the hospitality of the Institut 
fur Stochastik at the University of Freiburg and of the Isaac Newton 
Institute where part of this work was done. R.F. would like to thank 
the Institute for Advanced Study for its hospitality and support. 

2. Preliminaries 

When we generalize a problem in physics, usually additional degrees 
of freedom show up. In the case of the object one would like to call 
"general SLE" one faces similar questions. 

First, as we increase the connectivity of the domain, or look at ar- 
bitrary bordered surfaces, the possible simple paths of interest can do 
more, such that the notion of "radial" and "chordal" are not the only 
natural ones. 

Let us illustrate this with the example of the annulus. A simple 
curve, starting at the outer boundary could either terminate at an in- 
terior point, connect the same boundary component or end at the inner 
boundary. Therefore we have to make a choice at the very beginning 
by restricting the possible behavior of the curves. 

However, there is a second component. As we are ultimately in- 
terested in random paths on surfaces, we have to define a probability 
measure on them. Since this is not given by geometric considerations, 
additional data is needed. 

The natural point of view, at least from the standpoint of physics, is 
to see the curves as some characteristic manifestation of fields, defined 
on a (bordered) surface, which in the path-integral formulation, would 
be random, or quantum fields. A familiar geometric example of curves 
obtained via fields, are geodesies which are "manifestations" of the 
symmetric non-degenerate two-tensor field, i.e. the metric. The field 
perspective naturally connects with CFT, at least if we are interested 
in conformally covariant properties, and opens up the door for the 
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construction of a measure. But to obtain from a field theory the desired 
class (in the geometric sense) of random curves, we have to choose 
boundary conditions. For a continuum theory the natural ones are 
either Neumann or Dirichlet and combinations of the two. 

Again, in the familiar picture of statistical mechanics models on 
planer domains of some connectivity, such (random) curves would cor- 
respond to domain walls, connecting points where the boundary con- 
ditions change discontinuously. This is a "static" way of constructing 
measures on paths on surfaces, i.e. via the path-integral formalism 
(partition function). 

However, by the use of the Loewner mapping we can transfer the 
problem from the static point of view to a dynamic one by character- 
izing the arising curves via their driving function. 

Such a fundamental approach was first introduced by Oded Schramm 
in JH] in order to describe the scaling limit of two specific models, 
namely the loop-erased random walk and the uniform spanning tree. 

The natural assumptions for the above implied Brownian motion as 
driving function. 

But a closer look reveals, that the most natural driving function for 
models arising from say statistical mechanics, is to assume a Markov 
process, which allows for an additional drift. 

The crucial and new insight is, that the link between the set of 
boundary conditions, which generate the desired class of random curves, 
and the corresponding Markov process is given by representation the- 
oretic arguments. In the language of CFT, at the position where the 
boundary conditions change, we can insert boundary condition chang- 
ing operators, i.e. local fields \h) of some real weight h. 

As Hadamard's principle of boundary variation has a natural de- 
scription in terms of operators, which one can show to satisfy the com- 
mutation relations of the Virasoro algebra, the singular deformation of 
the boundary is given by the operators L_2 and L_i which correspond 
to "cutting" the surface and "moving" the point, where the boundary 
conditions change, i.e. the "marked point" . 

The condition on a vanishing drift as well as the way to calculate it 
at all, is exemplified in the following expression 



As we already mentioned in the introduction, another novel feature 
of the general Loewner evolution is the fact, that the driving function 
now lives on a higher- dimensional parameter space, i.e. a moduli space. 



( 



|l 2 _ x - 2L_ 2 ) ^ = . 
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Finally we remark, that for a multiply connected domain, we still 
could encode the random curves by a random driving function on a 
one-dimensional space, but then the resulting stochastic process would 
not be Markov anymore, i.e. the path would have to remember its 
complete past, and the stochastic process would be given by a Girsanov 
transform. 



3. Standard domains and canonical mapping 

Denote by D a domain in the complex z-plane bounded by a finite 
number of proper continua Cj,j = 1, . . . ,n, and let w be a point in 
D. If D is simply connected, then there is a unique conformal map 
= <&(z,w) from D onto the unit disk such that 

, . d$(z, w) , 

(1) $(w,w) = Q, Y z ' \ z=w >0, 

If n > 1, then there is a unique conformal map from D onto the unit 
disk with n — 1 disjoint concentric circular slits which maps C n to <9B 
and satisfies ((!]). We call the unit disk with a finite number of concentric 
circular slits a standard domain and the normalized mapping $(z, w) 
the canonical mapping. We now recall the construction of the canonical 
mapping. 

For z,w G D denote Gd(z,w) = G(z,w) the Green function of 
D with pole at w. G(z,w) is a harmonic function of z throughout 
D except at w, where G(z,w) + In \z — w\ is harmonic. Further, if z 
converges to a boundary continuum Cj, then G(z, w) converges to zero. 
These properties determine G. 

If D is simply connected, denote H (z, w) a harmonic conjugate of 
G(z,w) with respect to z. z H(z,w) is multiple-valued (as z de- 
scribes a small circle around w, H(z, w) changes by 27r) and contains an 
arbitrary constant. If we choose the constant such that for one branch, 
and x real, lim^o H (x + w,w) = 0, then the canonical mapping is 
given by 

(2) w) = exp(-G(z, w) - iH(z, w)). 

If n > 1, then G has also periods with respect to circuits around Cj, 
j = l,...,n- 1, 

(3) 27riu>j(z) = —i f ' ^ ds z . 

. r on,,, 



Here d/dn denotes the derivative in the direction of the outward point- 
ing normal and ds denotes arc- length measure. For the purpose of this 
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definition we have assumed that Cj is piecewise smooth. By the Rie- 
mann mapping theorem and the conformal invariance of the Green 
function this represents no loss of generality, uij(z) is harmonic in D 
with boundary values 1 on Cj, on Cfc, fc ^ j, and is called the har- 
monic measure of Cj in z with respect to D. It is also the probability 
that planar Brownian motion started in z exits D through Cj. Denote 
F(z, w) an analytic function in z with real part G(z, w). Then Uj(z) is 
the real part of the analytic function 

(4) R (z) = -- I 8F{Z ' W) ds 

Rj{z) is regular in D and possesses periods with respect to circuits 
around Ck given by 

(5) 2<KiP kj = i ! dSz . 

Jc k dn z 

It is well known that the period matrix P = [-P/y]/! is symmetric 
and positive-definite. Define the function 

(6) F(z,w) + R(z) t P- 1 u:(w), 
where 

~R(z) T = (Ri(z), . . .,R n _ 1 (z)), lo(w) t = (ui(w), . . 

This function has vanishing periods about Cj,j = 1, ... ,n — 1, and 
period — 2tt about C n . It is now easy to see that, after adding an 
appropriate imaginary constant, the canonical mapping for D is given 

by 

(7) $0, w) = exp(-F(z, w) - R(zfp- l u>(w)), 

see [in] for details. 

We can use $(z, w) to produce mappings onto other families of stan- 
dard domains, see 0. For example, if w — u + iv , then 

z i — > — — In w) 
ou 

maps D onto the whole plane slit along n segments parallel to the 
imaginary axis. The point z = w corresponds in this map to the point 
at infinity. Similarly, 

z i — ► — — In $(2;, w) 
ov 

maps D conformally onto the plane slit along segments parallel to the 
real axis. By letting w approach one of the boundary components, say 
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C n , we find that 



t/ \ d . . dF(z,w) ^ . T .,0u>(u>) 
u>) := - — In $(z, w) = — - - R(zY P 1 



dn w dn w dn 

is a conformal map from D onto the right half-plane 3?(C) > slit 
along n — 1 segments parallel to the imaginary axis. If the boundary 
components Cj are Jordan arcs, then \1/ extends continuously to the 
boundary. Note that if D is the unit disk and w = e llfi , then \1/ is given 
by 

+ z 

z i— > — . 

- z 

4. Domain constant and conformal radius 
We define the domain constant do{w) by 
(9) do(w) = — lim (G(z, w) + In \z — w\) , 

z— *w 

and the conformal radius rjj(w) by 



(10) r D (w) = In 



dz 



If D is simply connected then we find from (j2J) that dn{w) = td(w). If 
n > 1, then it follows from (J7J) that 

(11) rr>{w) = — lim (7(2, if) + In |z — w|) , 
where 

(12) 7(3, tu) := ^(ln$(^w)) = G(z,w) +w(w) T P" 1 w(z). 

We note that if D is a standard domain, then 7(2,0) = — In |z| and 
thus 

(13) G(z,0) = -\n\z\-u(0) T P- 1 u(z). 
It follows from (JUJ), ©, and JEJ that 

(14) d D {w) = r D {w) + u{w) Tn P~ l u{w). 



We call a closed and simply connected set A a hull in D if An D = 
A and B)\A has the same connectivity as D. Suppose that E is a 
standard domain and consider two disjoint hulls in E, say A and B, 
that do not contain 0. Denote *&a(z) = ^a(z, 0) the canonical mapping 
from E\A onto the standard domain E* := <&a(E\A) that fixes zero. 
Let D = E\B, D* = E*\$ A (B), and set T = dD, T* = dD*. We 
want to compare the "domain constant increments" do* (0) — do (0) and 
o?e*(0) — <i E (0). Assume that B is connected, intersects the unit circle 
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<9D, and has a piecewise smooth boundary. Since G E * (z, w)—G E * (z, w) 

is harmonic throughout D*, the Poisson formula gives 

(15) 

G D *{z,w)-G E *{z,w) = J {G D * (77, w) -G E * (r), w)) — |^ ^ ds v , 

where d/dn denotes the derivative along the outward pointing normal 
and ds integration relative to arc length. Note that we integrate along 
both sides of the slits. Since G* D (rj, w) = for r\ € T*, and G E *(r), w) = 
for rj G T*\$ A (B), we get fromjTSJ) and © 

(16) tMO) - MO) = ^ / G^faO) ^ fa0) rf v 
Similarly, 

(17) d D (0) - d E (0) = ± [ G E ( V , o) dG ^ V ' 0) ds v . 

^ JdBnE U n r) 

Let now {B e , e > 0} be a family of hulls as B above, such that B e 2> B £ > 
if e > e', and diam(B e ) = 0(e). Let £ = Q e>0 B e . We note that 
necessarily £ e In the following we suppress the subscript e when 
there is no risk of confusion. For rj e dB £ (1 E we have by Hadamard's 
formula 

dG D ( V ,0) dG E (£,0) 



dn v dn v 



+ o(l) 



dGB.($ A (ri),0) _ dGE-{*A(Z),0) | 
and by Taylor's formula 

^,0)=(-^M + O(e ))|,- { |, 

(19) G £ .(*^(.,),0) = (- aGE " ( ^ 4(0 '° ) + 0(e)) - <M0|, 
and also 

(20) / |t7-M0l <fc„ = |*k(£)| 2 / l^-el ^ + o(e 2 ). 

J<S> A (dBDE) JdBDE 

Thus, from (fTKl). (fTTjl . lfI8|l . IfEty . and lj2 ty . we get finally 

f2n rf^(0)-^(0) _/ dnv \ | $ w n |2 
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If the hulls B € are not smooth, then, given any 5 > 0, we may 
approximate them by smooth hulls B t ^, such that 

d Ds (o)-d E M <Mo)-<k,.(o) 

- cb(0)-d E (0) 'tfc.(0)-d B .(0) - + • 

Hence equation (J2TJ) applies also in the non-smooth case. 

Next, we wish to compare the "domain constant increment" with the 
"conformal radius increment." By (fT^j) 

(22) d D (0)-d E (0) = cj d (0) t P^lj d (0) - a^(0) r P^(0) 
1 ] r D (0)-r E (P) r D (0)-r E (P) 

We consider the domain D as a variation of the domain E, i.e. we con- 
tinue to consider the family {B e , e > 0} from above. From Hadamard's 
formula, [THl (45)] and [T3J Chap. I, Sec. 11], we have 

(23) Ifr/'W-'-sPQ-i 

e->0 £ 

(24) lim = 

e-^o e on orii 

(25 , lim P g -P E= g^).^«) T 

e^o e an an 

By (PU), it follows that 

(26) hm = — 1 W£;(0) P E — 5 . 

e->o e an on on 

We also have 

-1 r»— 1 a. a..!<-\T 



(27) lim P ^ ~ P * - -P-H) . W p-i 

,™ e _ ^ «9n dn E 



Combining and© we get 

fa rf g (0)-M0)^/ 1+ (ofp ,8^(1) 
r D (0) -r s (0) V 9" 

(28) /ae^or 

\ an 

From fl2TJ) and (j2HI) it now follows that 

(29) Um T "®- r ~'® = W A ®f 
e-o r D (0) - r E (0) 
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5. Radial Komatu-Loewner equation 

In this section we derive what we call the radial Komatu-Loewner 
equation. It describes the evolution of slit mappings when the slit 
grows from a point on the boundary of a multiply connected domain 
to a point in the interior. It shows in particular, that each growing 
slit induces a continuous motion on one boundary component of the 
domain, namely the image of the tip of the slit under the canonical 
mapping. The radial Loewner equation, of course, is the special case 
when the domain is simply connected. 

In |Zj , Komatu derives a differential equation that is satisfied by the 
canonical mappings for the growing slit in the case where the slit grows 
from one boundary component to another boundary component. The 
radial Komatu-Loewner equation is similar to the equation Komatu 
derives in [7] and our proof proceeds along the lines of the proof he 
gave. However, the equation given in |7] involves the derivatives of 
certain moduli with respect to the slit-parameter and it is only stated 
that these derivatives exist, but their explicit form is not given. Since 
we wish to solve the radial Komatu-Loewner equation for a given input, 
we need the explicit form of the equation. Using the fact that the real 
part of a single-valued analytic function is orthogonal to the real parts 
of Abelian differentials of the first kind, we derive the explicit equation 
directly. 

For the purposes we have in mind — generating random slits — the 
main difference between the simply and the multiply connected case 
is, that in the multiply connected case a growing slit corresponds to a 
motion on the boundary coupled to a motion on the moduli space. In 
the simply connected case the moduli space reduces to a single point. 

Let D be a standard domain, and 7 : [0, £ 7 ] — ► D a Jordan arc such 
that 7(0) G S 1 , and 7(0, t 7 ] C Z)\{0}. Let g t be the canonical mapping 
from D\7[0, t] which leaves fixed, and denote D t the standard do- 
main g t (D\y]f), t]). By ()23|). the map t 1— *• 5^(0) is a strictly increasing 
continuous function, see jTH] for details, and we may assume that the 
parameter t is chosen so that <^(0) = e*. We wish to find a differential 
equation for the family {g t : t G [0,£ 7 ]}. 

Denote Cj(t),j = l,...,n, the boundary components of D t . We 
always have C n {t) = S 1 . For j = 1, . . . , n — 1, let rrij(t) be the radial 
distance of the circular slit Cj(t) from the origin. Denote £(£) the 
starting point on S 1 of the Jordan arc gt(l[t,t 7 ]) in D t , i.e. gt(lt)- For 
< t* < t < t 7 , set 



9t,v = 9t* °9t X - 



ON RADIAL SLE IN MULTIPLY CONNECTED DOMAINS 



11 



Then g t> t* is a conformal map from D t onto D t *\g t * (7ft*, t]). The point 
= gt*{lt*) corresponds to two prime ends in D t *\g t *{^[t* , £]). De- 
note exp(z/? (t, t*)) and exp^/^t, £*)), with (3 (t,t*) < (5x{t : t*), the 
pre-images of these prime ends under g t>t * , i.e. 

g t>t *(exp(i(3 (t,t*))) = g t , t .(exp(i/3i(t, t*))) = gt*(lt*)- 

Then, if \z\ = 1 and ft (*,**) <argz< /3 Q (t,t*) + 2tt, 

\9tA*)\ = 1 - 



Consider the function 



In 



9t,t* (z) 



z 



Since g t>t * (0) 
(30) 



0, it is regular throughout D t . In particular, 

9t,t*(z) 



lim 

2^0 



j*-t 



Furthermore, as z describes a small circle around 0, In z changes to 
In z + 2ni and In g t>t * (z) changes to In g t>t * (z) + 2ni. Thus ln(g t)t * (z) / z) 
is regular and one-valued throughout D t . As In \gt,f(z)/z\ is regular 
and harmonic in D t , Poisson's formula gives 



(31) 



In 



9t,t*( z ) 



1 

2ix 



— / In 



dD t 



gtAQ 



C 



dG(C,z;t) 
dni 



ds, 



where G((, z; t) is the Green function for D t with pole at z. We write u 
for the harmonic function of z on the left-hand side of (|31|) and denote 
by v its harmonic conjugate. Since u is the real part of a one- valued 
analytic function, the period Pj(v) of v with respect to a circuit about 
Cj(t) has to vanish for each j = 1, . . . , n — 1. Since LUj(z; t) = 5j t k for 
z G Ck(t), we have 



= ^0;) 



dv 



Cj(t) 



c 3 {t) 



dv 

OS 



du , 

^~ ds 
on 



(32) 



, . du , 

LOj^t)— ds 

dD t dn 



u 



dD t 



Cj(t) 

dn 



ds. 



Combining 

9t,t* (z) 



and (|3*T|). we find 



In 

(33) 



1 

'2^ 



In 



8D t 



9t,AQ 



\ dni dn 



ds. 
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It follows from Section El that 



dG((,z;t) 



-iSw(C;*) 



dni ' dn 

has a single- valued harmonic conjugate and so 

9t,t*(z) 



In 



2tt 



/ In 



dD t 



9t,t*{0 



\ on\ On 



ds 



(34) 



+ ic, 



where z i— > F((, z; t) is a multiple-valued analytic function with real 
part G((,z;t), and where c is a real constant. Note that by (@J) and 
(JSJ), for every j = 1, . . . , n — 1, 



2tt 7 Cj(t) V dni 



dn 



ds 



(35) 



R j (z;t)-R(z;t) T p-\P t )., j = 



Thus, since In \gt,t*(()/(\ * s constant on each Cj(t), j = 1, . . . ,n — 1, 
and identically zero on {\z\ = l,(3i(t,t*) < aigz < fio(t,t*) + 2n}, it 
follows from (J53j) that 

9t,t*(z) 



(36) 



z 
1 

2^ 



ft (*,**) 



hi 



Po(t,t*) 



9tA0 



\ on\ on 



ds 



+ ic. 



We now show that c = 0. To that end, note first that by Cauchy's 
integral formula, 



(37) 



t* -t 



1 



In 



dD t 



9t,t*(0 \ d( 

c J c 



In particular, the right-hand side of (|3T|) is real. Since all boundary 
components are concentric circular slits, d(/( is purely imaginary along 
dD t , i.e. 

d( 



C 



i darg(, ( G dD t 
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Hence 



t*-t 



2tt 



In 

dD t 

(38) = - / In | ^ ( C *) \dip+-'£ ^ ! ^77V d ar ^' 



<M0 



<iarg£ 



In 



Since the two "sides" of Cj(t) make opposite contributions 



dargC = 0, j = l,...,n- 1, 



C,-(t) 



and we finally get 



(39) 



t* -t = — 
2tt 



A (*>**) 



lnk,,(e^)| dp. 



A>(t,**) 



Note next that, by (J7J), 



F(C, 0; t) + R(C; tfP^"^ t) = - In C 

Further, since we are free to add a constant to the harmonic conjugate, 
we may assume that 3(i?(0;t)) = 0, i.e that u>(0;i) = i?(0;£). Also, 
since u>(£;t) is constant along the boundary, 



dn 

Thus, using P = P T , we get 



91nC 

dn 

drii 



ds 



0. 



(40) 



+ P(0;tfp-^ 



Hence, if we evaluate (J3fij) at z = and use (J3Uj) on the left, and (J3HJ), 
(j4*Uj) on the right, then it follows that c = 0. 
Letting z = gt(w) we now get 



In 



9t*{w) 



9t(w) 



2^ 



(41) 



ln\g t , t 4e^ 
dF(e ilp , z; 1) 



x 



dni 
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We now wish to let t* / t in (JUJ). Note first that, for cp G [0,27r], 
<p i— > In | (e* 1 ^) | is continuous and non-positive and that also 

A , N dF(e^,z:t) n , , T „ 1 9w(e i *';*) 

is continuous. Thus it follows from the mean-value theorem of integra- 
tion that 

i rlh (<>**) 

mk t »(e^)U(^) dtp 



2n(t*-t) J mr) 

$t(A(<p')) + i$S(A(<f/')) f Mt ' n 



/ ln|^(e^)| d 9 

Jdnft.t*) 



(42) = 8^(^))+i9(A(^)), 
for some (p',(p" G [/3 (t, t*), t*)\. Hence 

lim hig t *(w) - \ng t (w) 

(43) 7 = jhk^ll _ jK^tf^MfW;*), 

ani on 
By the same argument we may let t \ t*. On the right-hand side 
above we then only need to change t to t* and introduce an overall 
minus sign. Thus we have established the following 

Theorem 5.1 (Radial Komatu-Loewner equation). If 7 is a Jordan 
arc in a standard domain D starting on S 1 with the parametrization 
from above, and if g t is the canonical map for _D\7[0,t], then, using the 
notation from above, the family {g t : t G [0,i 7 ]} satisfies the equation 

(44) ft]n* W = .Bl*).^^^, 

a^i an 
with initial condition go(z) = z. 

6. Motion of moduli 

The right-hand side of the radial Komatu-Loewner equation, at time 
t, involves the Green function of the domain D t , and also various 
functions derived from the Green function. Consequently, it does not 
make sense to ask for the solution of (|44jl for a given continuous curve 
t I— > £(£), since the vector-field on the right-hand side of (}4*4*|) is not 
specified by giving that information alone. To specify the Green func- 
tion of D t we also need the moduli of the domain D t . We will now 
consider what the appropriate moduli space is for our purposes and 
find a system of equations these moduli satisfy. Once this system is 
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found, we can solve it for a given input t i— > £(t), and then, in a sec- 
ond step, solve the radial Komatu-Loewner equation using £ and the 
moduli. 

The geometric description of D t requires 3n — 3 real parameters, 
three for each (interior) slit, given, for example, by the radial distances 
of the slits to the origin, i.e rrij(t), j — 1, . . . , n — 1, and the angles 



determining the endpoints of the slit Cj(t), j = 1 . . . ,n — 1. On the 
other hand, it is well known that two n-connected domains with non- 
degenerate boundary continua are conformally equivalent if 3n — 6 real 
parameters agree for n > 2. If n = 2 then there is only one real 
parameter describing the conformal class, and if n = 1, then all such 
domains are conformally equivalent. In our context, we only allow 
conformal maps for which a certain interior point has a prescribed 
image and whose derivative at that point is positive — the canonical 
maps from Sectional This corresponds to considering domains with one 
marked interior point and one marked boundary point. Indeed, given a 
domain D with n non-degenerate boundary continua, one interior point 
z and one boundary point (or, more generally, prime end) (, there is a 
unique conformal map from D onto a standard domain such that z is 
mapped to 0, and ( to 1. 

For the slits we wish to grow the marked points are the beginning 
(t = 0) and end points (t = oo). It is now easy to see that the 
moduli space of n-connected planar domains with one interior and 
one boundary point marked is 3n — 3 dimensional for all n > 1. 
We will take m(t) = (m^t), . . . , m n _i(t)), 6{t) = (9t(t), . . . , n _i(*)), 
and 6' {t) = (0[(t), . . . ,6'^_ 1 (t)) as the moduli of the domain D t and 
write M(t) := (m(t) , 6 (t) , & (t)) . For a standard domain the marked 
points are and 1. For a point M in the moduli space we denote by 
D = D(M) the corresponding standard domain, while for an arbitrary 
n-connected domain D we write M = M(D) for the corresponding 
point in the moduli space. 



Then, by (jHJ), maps D t conformally onto a slit right half-plane. By 
(HUJ), ^ t is the unique such map with 



< J = l,-.., 



n — 1, 



Set 



(45) 




(46) 
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We may write the radial Komatu-Loewner equation (|44jl as 

d t \ng t (z) = %(g t (zU{t)). 

By boundary correspondence, if z G Cj, then gt{z) G Cj(t) and 

»(ln<7 t (z)) = lnm i (t). 

Thus, by considering the real part of the radial Komatu-Loewner equa- 
tion, 

(47) ^lnm i (t) = K(* t ^),e(t))). 

Further, if 

are the endpoints of the slit Cj(t), then 

Zj {t) = ^K-(O)e^W), *>(f) = ^K-(O)e^W), 

where 9j(0) < T)j{t),r)j(t) < 6j(0). Indeed, the pre-images of the tips 

of Cj(t), that is 7n J -(0)e M7 -'( t ) and m^C^e 1 * 7 ^, are the solutions to the 
equation 

on the set of prime-ends corresponding to Cj\{zj(0), z'j(0)}. A tip of 
Cj(t) cannot be the image of a tip of Cj because then the analytic 
function dgt/dz would not have the required number of zeroes, 2n — 2. 

Lemma 6.1 (Motion of moduli). The moduli 

M(t) = (m(t),6(t),0'(t)) 

satisfy the system of equations 

d t ln mj (t) = - ( P; 1 



dn 

d t j (t) = Q(% ( mj (t)e i0 M,at))), 
(48) ^(t) = 9f(* t (m i (t)e^W 

for j = l,...,n- 1. 

Proof. We note that dg t /dz and d 2 g t /(dz) 2 are analytic functions that 
extend analytically to the prime-ends corresponding to Ci,...,C„_i 
with the endpoints of the slits removed. By the implicit function the- 
orem, 
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is differentiable with derivative 



DER t 



^( mi (0)e*w«; 



(dz 

By counting zeroes we find that 



d 2 9t 



WO)^**) ^ 



(dz) 2 

and so DER t is finite. Hence 

^■(t) = 9t(9f(^(m i (0)e^W)) 

= 3 (% (m 3 (t)e l6 ^\i(t))) + ^ x DER t ) 

(49) =9f(* t (m i (*)e^W,e(*)))- 

In a similar way we obtain the derivative of ^-(t). It remains to check 
that ()47|) agrees with the first equation in (|48p . To this end we note 
that 

R(»,(,,0) = -^|f^ - w( ,;^p,-'^. 

From the boundary behavior of the Green function and the harmonic 
measures, it follows that for z G Cj(t) 

dG((,z;t) 

= 0, and u k (z) = b jk . 

dni 

The lemma follows. □ 

We now have our main existence statement. 

Theorem 6.1. Given a continuous function t G [0, oo) i— > £(t) G S 1 
and the moduli M of a standard domain D, there exists a unique so- 
lution M(t) to the system (|48|) on an interval [0, t^) with iVf(O) = M, 
and where t% is characterized by 

t% = inf{r : limm^t) = 1 for some j G {1, . . . ,n — 1}}. 

t/r 

Further, if D t is the standard domain determined by M(t), and if 
*& t (z,C) is the holomorphic vector field associated to D t by (145)) . then, 
for any z G D, the equation 

d t \ngf(z) = %{gf(z)^(t)), g°(z) = z, 



has a unique solution on [0, t z ), where 



t z = sup{t < % : inf \g»(z) - £(s)| > 0}. 
se[o,t] 
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Finally, for t < t% set K t = {z G D : t z < t} . Then g® is the canonical 
conformal map from D\K t onto D t which fixes zero and has positive 
derivative there. 

Proof. For the existence of the solution to the moduli equations ([48)1 on 
[0, we will show that the vector field in ()48|) is Lipschitz as a function 
of M, with a Lipschitz constant that only depends on distance to £(t) 
of the slit (or slits) nearest to £(£). 

Let M and M* be two points in moduli space with corresponding 
standard domains D and D*, such that 

Im^-m*!, 1^-0*1, |0J-0?| <e. 

We assume that e is so small that 

Cj H C* k = 0, whenever j ^ k. 

Denote Zj, z* the endpoints of the slit Cj and z*, z'* the corresponding 
endpoints of C*. Denote ^ the vector field for D and ty* the vector 
field for D*. Then we need to show that 

(50) t!*{z*) - *( Zj ), V*(z'*) - = 0(e), ./ 1 // 1. 

We begin by applying a particular interior variation as in Garabe- 
dian's proof of Hadamard's variational formula, [3]. 

Denote xi,X2 real coordinates for D. It is easy to see that we can 
map D one-to-one onto D* by a transformation 

(51) Xj^Xj + eSj, J = 1, 2, 

which sends the endpoints Zj, z'j of the slits in D to the endpoints z*, z'*, 
respectively, of the slits in D*, and where 

Sj = Sj(xi,x 2 ), J = 1, 2, 

is a pair of smooth functions in some neighborhood of the closure of 
D. In complex notation, 

z* = z + eF(z, z), 

where 

F(z, z) = Si(x 1 , x 2 ) + iS 2 (xi 1 x 2 ). 

Similarly, we put C* = C + e -^(C) 0? w * — w + eF(w, w). Denote 

g(z,(;e)=G*(z*,w*) 

the transformed Green function. Here G* is the Green function of 
D*. Then it is straightforward to show that g satisfies the equation 
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L e [g] = 0, where L e represents the self-adjoint differential operator 
defined by 

2 



k,l=l 

with coefficients 

-1 2 



= / gCf^g) \ y- ^ ^ 



d(x 1 ,x 2 )J ^dxkdxi' 

Note that A n A 2 2 — A 12 A 2X = I and that 

(52) A fei = S ki + z<t>ki (xi ,x 2 ,e), 

where <pki is smooth for x±,x 2 in a neighborhood of the closure of D 
and e in a neighborhood of zero. Let a,jk = djk(z) be the inverse matrix 
of A, and denote T the quadratic form 

2 

j,k=l 

where £i and £2 stand for the real and imaginary parts of (. Finally, 
denote a = a(z, () a fixed smooth function of the two points z and ( 
in D, which fulfills the boundary condition a(z,Q = when either z 
or ( lies on dD, but has the value 

when z coincides with ( inside D. Then 



defines a parametrix for L t . Direct computation shows that 

L t [P e ]-L [P ] = e<j ) (zX,e), 

(53) L e [G - P ] -L [G- P ] = ei>{z, C, e), 

where and ip are smooth functions in z,(, and e, except for a 1/r- 
singularity on the diagonal (here r = \z — (\). Using Green's identity 
and the boundary conditions, 

g((, w; e) = G(w, () - P (w, () + P £ (w, () 

(54) + / g(z, w; e)L e [G - P + P e ] dx 1 dx 2 , 

JD 



20 



ROBERT O. BAUER AND ROLAND M. FRIEDRICH 



see We note that the integral exists in the sense of Lebesgue as the 
1/r singularity is integrable. It then follows from (J52J) and that 

g((,w;e)-G(wX) = 0(e), 

uniformly in (, w in the closure of D if \(— w\ is bounded away from zero. 
If we take the normal derivative of both sides of (|54|) with respect to 
w G dD, then the identity continues to hold as the l/\z—w\ -singularity 
that now appears in the integrand is still integrable. We find 

- g ((, W ;e)- 1 —G(w,() = 0(e), 



dn w <9n u 

uniformly in ( G D if \( — w\ is bounded away from zero. If we take a 
further derivative in ()54j) . 

d_ _ 1 ( d_ _ . d 

then the integrand is no longer integrable in the sense of Lebesgue and 
has to be understood in the sense of a Cauchy principle value. Thus, 
(IH4l . together with (J52J) and (|53J) still imply 

< 55) iS^j^-iS^ G(a ^= o ^ 

uniformly in £ G 5 if |£ — w\ is bounded away from zero. We now 
note that all domain functions that are used in the construction of the 
vectorfield ^ can be obtained via integration from d 2 / (d(dn w )G(w , (). 
Thus we will be done if we can show that 

d 2 d 2 

-G*(e,w*)-——G(w,() = 0(e), 



d(*dn w * d(dn v 

uniformly in £ G D if \Q — w\ is bounded away from zero. But this 
follows from and 

d 2 d 2 

-G*(e,w*)-——g((,w,6) = 0(e), 



d(*dn w * d(dn h 

the latter being a consequence of (joTj) . 

The second part of the theorem now follows from general results 
about ordinary differential equations, exactly as in the simply con- 
nected case. □ 
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7. Radial SLE in multiply connected domains 

The purpose of this paper is 1) to give a "natural" construction 
of conformally invariant measures on "simple curves" in multiply con- 
nected domains, and 2) to study some of the properties of these random 
curves. We will now motivate, using informal arguments, our partic- 
ular construction of conformally invariant measures on simple curves. 
The arguments lead to a small class of processes which contains radial 
SLE K in multiply connected domains. 

For a domain D with n non-degenerate boundary continua, a bound- 
ary point z G dD, and an interior point w G D, let W(D, z, w ) be the 
set of Jordan arcs in D with endpoints z and w. Denote {C^ z w }d,z,w 
a family of probability measures on Jordan arcs in the complex plane 
such that 

C% z jW(D,z,w)) = l, 

and where M = M(D). Such families arise, or are conjectured to 
arise, as distributions of interfaces of statistical mechanical systems at 
criticality. Based on these models, e.g. percolation, one expects that 
the distributions describing the interfaces in different domains with 
different marked points are related by a Markovian-type property and 
conformal invariance. Denote 7 a random Jordan arc with law C^ zw . 
The Markovian-type property says that if 7' is a sub-arc of 7 which 
has z as one endpoint and whose other endpoint we denote by z', and 
if M' = M{D\ r y') ) then the conditional law of 7 given 7' is 

(56) law( 7 |7') = £^' 7V> . 

Conformal invariance means that if / : D — > D' is conformal, z' = f(z), 
w' = f(w), then 

(57) ^D 1 ,z' ,w' = f*£i?,z,w 

If (j5T)l holds, then to understand the family {C^ z w } it is enough to 
consider standard domains D, take w = 0, z = 1, and, by the identifi- 
cation of standard domains with their moduli, we may write 

rM _ pM 

In this case there is a natural parametrization of the Jordan arcs we 
consider. Let 

s G [0, 00) 7(5) G D 
be a Jordan arc in a standard domain D such that 

7(0) G S\ 7(0,00) C D\{0}, and lim 7(f) = 0. 
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Denote M = M(D) the point in the moduli space corresponding to D 
and let gf 1 be the canonical mapping from .D\7[0,t] onto a standard 
domain D t := g™(D\y[0,i\). Then 

W(0)<(rfY(0), 

whenever t* < t, and lim^oo (g^ 1 )' (0) = oo. Thus we may and always 
will assume that 7 is parametrized by the conformal radius of S 1 U 
7(0, •] in with respect to dD, i.e. so that (g^ 4 ) (0) = e*. This 
parametrization is natural in the following sense. If t > 0, M(t) = 
M(D t ), and 7 is the curve defined by 

se[O,oo)H+7( s )=0 t M ( T (t + s )), 

then the canonical mapping gi^'^ from D t \7[0, s] is given by 

9 f (i) =^°U M )" 
and so #s ^(.D t \7[0, s]) = -D t+S . In particular 

b. M(t) ) / (0) = e^e-* = e ' I 

i.e. 7 is also parametrized by conformal radius. 

Let now {g^ 1 : s > 0} be the random family of canonical maps 
corresponding to the random Jordan arcs {7[0, s] : s > 0} in a standard 
domain D, and denote 

C M = law({^ M : s > 0}). 

Then, applying first the Markovian-type property and then conformal 
invariance, (JEHJ), (|H7j) . we find 

iaw({^ :s >o}i^) = (,n; 1 ^ w - 

Equivalently, 

(58) law ({«& o : a > 0}|<?f ) = law({^° = « > 0}). 

By the radial Komatu-Loewner equation, (|44j). for each t > 0, the <r- 
field generated by g t M is equal to a((9(r), M{r)) : r G [0, t]), where 
exp(i#(r)) = £(r), and 0(0) = 0. Similarly, it is easy to see that we can 
reconstruct g£* s o (g^)' 1 from {(0(t + r) - 9(t), M(t + r)) : r G [0,s]}. 
Thus (J55j) implies 

law({(0(* + s) -0(t),M(t + s)) : s > 0}\{(9(r), M(r)) : r G [0,t]}) 

(59) = law({(0(s), M(s)) : s > 0}), 

where M(s) = M(D t \j[0, s]), for a random Jordan arc 7 with law 
£Af(t)_ The equa iity (j^J j s precisely the statement that {(0(t), M(t)) : 
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t > 0} is a Markov process. We note that in the simply connected case 
(n = 1), (jSUjl reduces to 



law({0(* + s) - 6{t) : s > 0}\{6{r) : r G [0,t]}) = law({0(s) : s > 0}), 



from which it follows that 9 is a process with independent, and identi- 
cally distributed increments. From this, continuity, and the symmetry 
law(0) = law(— 9), Schramm derived in JH] that 9(t) = \^B t for a 
standard one-dimensional Brownian motion and a positive constant 
k. The continuity follows from the continuity of the Jordan arcs, and 
the symmetry is actually observed in various discrete models, such as 
loop-erased random walk. 

In the multiply connected case, we also have continuity. Let us now 
study the Ito differential for the diffusion (9(t), M(t)). By Lemma 
M(t) is a finite variation process and its differential is given by ([48)1 . 
On the other hand, we know from JS] that the qualitative properties 
of 7(t) would change with t if the martingale part of 9(t) has quadratic 
variation which is nonlinear in t. Thus d9(t) = ^JTidB t + "drift", and 
the only open question concerning the diffusion (9(t), M(t)) is the drift 
of 9(t). For a general drift, which may be a function of M, 9, and k, 
the resulting family {gjr*} is a random Loewner chain and we call the 
diffusion a Schiffer diffusion. In the case of percolation, the drift is 
easily identified. Consider a honeycomb lattice-approximation to our 
domain D. The beginning of the random simple curve 7 is an edge e 
of the lattice in D, which is horizontal and whose one endpoint is the 
point 1. The next step is either up or down (with slope 2ir/3). For 
percolation, each of these possibilities has probability 1/2. Denote p 
the endpoint in D of the slit made up of two edges, if the second step 
went up, and q the corresponding endpoint, if the second step went 
down. Denote g the canonical map from D\{e}. To derive the drift 
from this condition, we compare the images of the endpoints p and q. 
Recall that for the unit disk D 



F(z, w) = ln(l — zw) — ln(w — z). 



Hence 



dF(z, w) 



z + w 



z e s\ 



dn z 

For a standard domain D, let 



z — w 



k D (z,w) 



dF D (z,w) 
dn z 



z + w 



z — w 



ze S 



1 



and define 



ko{z) = lim kn{z,w), z E S 1 . 
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If the lattice size is small, say a single edge is length y/e, then g'(0) = 
e + o(e). Thus 

2 = 6 H 2 + ° (e) 

(60) = e(k D (l) + ^(l)^ 1 ^) + o(e). 

Hence, to model cluster-boundaries of percolation in a multiply con- 
nected domain D we make the ansatz 

do{t) = - 1 (k(m m)) + m(ty, ^)) r P4 ^ ( y (f)) ) * 

(61) + dBt, 

with f (t) = e ie(t) , and where M(t) satisfies (gHJ). 

For other discrete models the same reasoning as above would lead to 
different drifts. For example, for loop-erased random walk, the proba- 
bility of stepping up is not the same as stepping down, depending on 
the configuration of the concentric circular slits. However, the respec- 
tive probabilities can be calculated in terms of harmonic measure. This 
leads to a different Schiffer diffusion and we expect the resulting family 
{g®} to be closely related to the "harmonic random Loewner chains" 
studied by Zhan in his thesis, [T^J . We leave the question of which 
Schiffer diffusion corresponds to which discrete model to a forthcoming 
paper. In principle, representation-theoretic considerations as are done 
in conformal field theory should identify the relevant class of Schiffer 
diffusions. 



8. Locality 

In this section we show that the ansatz (jpTj) leads to random growing 
compacts satisfying the locality property if k = 6. Denote {gf, t > 
0} the solution of the radial Komatu-Loewner equation in a standard 
domain E starting at z = 1 for the Schiffer diffusion ([61]). Denote 
{K t ,t > 0} the associated growing compacts. Let A be a hull in E 
that contains neither zero nor z = 1. For the following calculations we 
restrict to the event {t < r}, where r := inf{t : K t n A ^ 0}. Let $a 
be the canonical mapping from E\A, g\ the canonical mapping from 
$a(E\(A U Kt)), and h t the canonical mapping from g t (E\(A U Kt)). 
Since the canonical mapping for E\(A U K t ) is unique, we have 

(62) h t og t = g* t o<S> A . 

Furthermore, up to a time change, the family {g t *} also satisfies a radial 
Komatu-Loewner equation beginning with the standard domain E* : = 
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$ A (E\A). In fact, it follows from flUJ) that 
(63) 

= \Kmm (MM + r*«; (pr 1 ^™ 

where u>^ = g^(z), and = h t (^(t)). The question we are interested 
in is whether (£*, iVf *) is a time change of (£, M). Since /ij = <j£ o <fe A o 
g t _1 , we have 

(64) dthtiz) = [d t g*] ^ A {9i\z))) + {g* o <$> A )\ g ;\z)){d t g;\z)), 

and we note that 
(65) 

Let t (^) = — i In h t (e lz ). Then d t 4> t (z) = — id t h t (e tz ) / h t (e lz ) , and 
#(z) = e i2 /4(e i2 )//i t (e i2 ). Thus (0),©, and (JEHJ) imply 

dtMz) ~ i dn 

+ ^W R > ( ^ W;t) T (pg) -i^(r(t);t) 
z an 

(66) ^ ( ' } * 5n 

Hence the stochastic differential 

d t <p t {z) dt + (f>' t (6(t)) d9{f) 

has martingale part <p' t (9(t))y/K dB t and its drift part consists of the 
three components 

mt) )> ^(rw,«*<M_ t . r(t);t ,V 



dn 

II : = mt)? ( k'(V(t); t) + R*(e*W; f (Kj-i ggfjjj 
z \ on 

(67) 

m : = «HW A. ({(i);t) + R(e , lfP -,M(M) A , 

z \ on J 
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When z — > 0(t), then part III converges to zero, and part II, together 
with the martingale part, converges to a time-change of (J61j) starting 
at E*. Finally, the limit of part / is by the definition of k(£;t) equal 
to 

(68) Um f-iw _ ?m) = - me) . 

z-^e \<f)(z) - (j){9) z -9 J 
Thus, by Ito's formula, 



i \ ' On 

(69) + H^.<%(0(t)) dt + <f>' t (6(t))Vt dB t , 

which is indeed a time-change of ()61|) if and only if k = 6. From f)6>3j) it 
follows immediately that the equations for M* are given by the same 
time change of the equations for M. 

Theorem 8.1 (Radial SLEg). The solution to the radial Komatu- 
Loewner equation based on the Schiffer diffusion (|61jl satisfies the lo- 
cality property if and only if k = 6. 
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